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Abstract 

We study the spherical collapse model in the presence of quintessence with negligible speed of sound. 
This case is particularly motivated for w < — 1 as it is required by stability. As pressure gradients 
are negligible, quintessence follows dark matter during the collapse. The spherical overdensity be- 
haves as a separate closed FLRW universe, so that its evolution can be studied exactly. We derive 
the critical overdensity for collapse and we use the extended Press-Schechter theory to study how 
the clustering of quintessence affects the dark matter mass function. The effect is dominated by 
the modification of the linear dark matter growth function. A larger effect occurs on the total mass 
function, which includes the quintessence overdensities. Indeed, here quintessence constitutes a third 
component of virialized objects, together with baryons and dark matter, and contributes to the total 
halo mass by a fraction ~ (1 + w)QQ/Q m . This gives a distinctive modification of the total mass 
function at low redshift. 

1 Introduction 

One of the most important open questions for cosmology is whether dark energy is a dynamical 
component of the universe or a cosmological constant. A plethora of experiments and future obser- 
vations are currently planned with the aim of improving our understanding of this question (see for 
instance pQ for a review). One of the most popular models of dynamical dark energy is quintessence, 
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where the acceleration of the universe is driven by a scalar field with negative pressure. Standard 
quintessence is described by a minimally coupled canonical scalar field [2J. In this case, scalar fluc- 
tuations propagate at the speed of light and sound waves maintain quintessence homogeneous on 
scales smaller than the horizon scale [3] . Quintessence clustering takes place only on scales of order 
the Hubble radius, so that its effect is strongly limited by cosmic variance. However, quintessence 
modifies the growth evolution of dark matter through its different expansion history. Thus, future 
galaxy catalogs and weak lensing surveys will have a great potential in constraining the dark energy 
properties, in particular its equation of state, through a detailed study of the evolution of dark 
matter. 

The study of quintessence perturbations - where by quintessence we indicate a dark energy 
sector described by a single scalar degree of freedom - has been recently revived in |4j. Here the 
most general theory of quintessence perturbations around a given background was derived using 
the tools developed in [5j [6], formulated in the context of an effective field theory. An important 
conclusion is that quintessence with an equation of state w < — 1 can be free from ghosts and 
gradient instabilities OH]. In this regime, stability can be guaranteed by the presence of higher 
derivative operators [3] with the requirement that the speed of sound of propagation is extremely 
small, |c s | < 10~ 15 . On cosmological scales, higher derivative operators are phenomenologically 
irrelevant and quintessence simply behaves as a perfect fluid with negative pressure but practically 
zero speed of sound [4j. Interestingly, this description applies also in the non-linear regime, i.e. when 
perturbations in the quintessence energy density become non-linear, as long as the effective theory 
remains valid. 

This study motivates the possibility that quintessence has a practically zero speed of sound. 
Apart from these theoretical considerations, the fact that the speed of sound of quintessence may 
vanish opens up new observational consequences. Indeed, the absence of quintessence pressure 
gradients allows instabilities to develop on all scales, also on scales where dark matter perturbations 
become non-linear. Thus, we expect quintessence to modify the growth history of dark matter not 
only through its different background evolution but also by actively participating to the structure 
formation mechanism, in the linear and non-linear regime, and by contributing to the total mass of 
virialized halos. 

In the linear regime, a series of articles have investigated the observational consequences of a 
clustering quintessence. In particular, they have studied the different impact of quintessence with 
c s = 1 or c s = on the cosmic microwave background [HI [10l [11] , galaxy redshift surveys [12] , 
large neutral hydrogen surveys [13], or on the cross-correlation of the integrated Sachs- Wolfe effect 
in the cosmic microwave background with large scale structures |14[ |15j. On non-linear scales, the 
dependence of the dark matter clustering on the equation of state of a homogeneous quintessence, 
i.e. with c s = 1, has been investigated using N-body simulations in a number of articles (see for 
instance [16J and references therein for a recent account). 
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A popular analytical approach to study non-linear clustering of dark matter without recurring to 
N-body simulations is the spherical collapse model [T7J. In this approach, one studies the collapse of 
a spherical overdensity and determines its critical overdensity for collapse as a function of redshift. 
Combining this information with the extended Press-Schechter theory [18|. [T9] one can provide a 
statistical model for the formation of structures which allows to predict the abundance of virialized 
objects as a function of their mass. Although it fails to match the details of N-body simulations, 
this simple model works surprisingly well and can give useful insigths into the physics of structure 
formation. Improved models accounting for the complexity of the collapse exist in the literature 
and offer a better fit to numerical simulations. For instance, it was shown in [20J that a significant 
improvement can be obtained by considering an ellipsoidal collapse model. See also |2H [22] for 
recent theoretical developments and new improvements in the excursion set theory. 

The spherical collapse can be generalized to include a cosmological constant (see for instance 
[23] ) and quintessence with c s = 1 [23] (see also [251 ES] for subsequent applications) . If quintessence 
propagates at the speed of light it does not cluster with dark matter but remains homogeneous. 
Indeed, pressure gradients contribute to maintain the same energy density of quintessence between 
the inner and outer part of the spherical overdensity. A study of the spherical collapse model with 
different quintessence potentials was performed in [27]. For a nice review on structure formation 
with homogeneous dark energy see also [28J. 

In this paper we study the spherical collapse model in the case of quintessence with zero speed 
of sound. This represents the natural counterpart of the opposite case c s = 1. Indeed, in both cases 
there are no characteristic length scales associated to the quintessence clustering and the spherical 
collapse remains independent of the size of the object. For this study we describe quintessence using 
the model developed in [2]. As explained, this description remains valid also when perturbations 
become non-linear and can thus be applied to the spherical collapse model. 

As the spherical collapse occurs on length scales much smaller than the Hubble radius, we will 
describe it using a convenient coordinate system where the effect of the Hubble expansion can be 
treated as a small perturbation to the flat spacetime. In these "local" coordinates the description 
of the spherical collapse becomes extremely simple and the well-known cases can be easily extended 
to quintessence with c s = 0. In this case pressure gradients are absent and quintessence follows 
dark matter during the collapse. Thus, in contrast with the non-clustering case c s = 1 where 
quintessence and dark matter are not comoving, for c s = the collapsing region is described by an 
exact Friedmann-Lemaitre-Robertson- Walker (FLRW) universe. Note that even though the energy 
density of quintessence develops inhomogeneities as long as the collapse proceeds, the pressure inside 
and outside the overdense region remains the same. Thus, as explained below, our model does not 



1 The characteristic length scale associated to the quintessence clustering is the sound horizon scale, i.e., 
L s = a J c s dt/a. As mentioned above, this vanishes for c s = so that clustering takes place on all scales. For 
c s = 1 we have L s = 2Hq , which is much larger than the scales associated to the spherical collapse. 
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give the same description of clustering quintessence as that proposed by [27] and studied, for instance, 

in [221 EH EB- 

We will see that, besides quantitative differences with respect to the c s = 1 case - a different 
threshold for collapse and a different dark matter growth function - the c s = case has a remarkable 
qualitatively new feature. Quintessence clusters together with dark matter and participates in the 
total mass of the virialized object, contributing to their gravitational potential. 

The plan of the paper is the following. In section [2] we describe quintessence models with c s = 0. 
In section [3] we study spherical collapse solutions first in known cases (dark matter only, ACDM and 
c s = 1 quintessence) and then in the case of a clustering dark energy, c s = 0. It turns out to be much 
simpler to describe these solutions in coordinates for which the metric is close to Minkowski around 
a point in space. The equation for the evolution of the spherical collapse are solved in section H] and 
the threshold for collapse is calculated in the various cases. This leads to the calculation of the dark 
matter mass function in section [5j In section [6] we study the accretion of quintessence to the dark 
matter haloes and its effect on the total mass function. The contribution of quintessence to the mass 
may be distinguished from the dark matter and baryon component in cluster measurements, as we 
briefly discuss in section [7J Conclusions and future directions are discussed in section [51 



2 The model: quintessence with c s = 

Let us consider a /c-essence field described by the action [32J [33J 

S = J d A x^g~ P(<P, X), X = -g^d^SA . (1) 



The evolution equation of <j) derived from this equation is 

1 



=0„(V=£ 2P tX d»ij>) = -Pj , (2) 



where Pj = dP/df. The energy- momentum tensor of this field can be derived using 

6, 



and can be written in the perfect fluid form as |34j 

Tpv = (pq + PQ)Ufj,u v + pQg^u , (4) 

once we identify 

PQ = 2P X X-P, PQ = P, u„ = -M (5) 

VA 



Let us initially neglect perturbations of the metric and assume a flat FLRW universe with metric 
ds 2 = —dt 2 + a 2 (t)dx 2 . The energy-momentum tensor of the field can be perturbed around a given 
background solution (f>{t) corresponding to a background energy density and pressure, 

PQ = 2P X X -P, Pq = P, (6) 
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where X = 4> 2 . To describe perturbations it is useful to write the scalar field as 011] 

= $(t + n(t,x)) , (7) 

where tt describes the difference between the uniform time and scalar field hyper surfaces^] Then, 
eq. ([5]) can be expanded linearly in ir using (p(t, x) = (f> + <p tt and X(t, x) = X + Xir + 2Xn. This 
yields, for the perturbations of the energy density, pressure and velocity, 

&PQ = h 71 + (PQ +PQ+ 4M 4 )TT , 5p Q = pQTT + (p Q + Pq)tT , Ui = -<9j7T , (8) 

where we have used eq. ^ and defined M = P jXX X 2 , where M has the dimension of a mass. 

To describe the evolution of perturbations we can expand the action (pQ) up to second order in ir 
as done in [3], 



S= d^xa^ 



Pit + 2P jX Xtt + (P jX X + 2P XX X 2 ) ir 2 - P X X^^- + \Pk 2 + 2 {P X X) 



a 2 2 

'(9) 

The second term proportional to n can be integrated by parts so that the part of the action linear 
in 7r can be written using eq. © as — [pq + 3H(pq + Pq)\ 7r, where H = a/a is the Hubble rate. 
This part cancels due to the background equation of motion. Furthermore, we can manipulate the 
last two terms of the action integrating by parts the last term, proportional to ittt, and making use 
of the background equation of motion, to rewrite them as 3HP jX X tt 2 . Finally, it is convenient to 
rewrite the coefficients left in this expansion in terms of the background energy density and pressure 
using eq. ©. This yields 



S = / d 4 x a 3 



~ (PQ +PQ + 4M 4 ) n 2 - ^(p Q+ p Q )^- + 1h{Pq +Pq)7t 2 



(10) 



The coefficients of this quadratic action are completely specified by the background quantities 
PQ +PQ an d M A . The latter is a function of time which we expect to vary slowly with a rate of order 
Hubble@ As shown in [6l H] , eq. ()10p is the most general action describing quintessence in absence 
of operators with higher-order spatial derivatives. Note that this action is even more general than 
the starting Lagrangian ([T|) as it can be generically derived using only symmetry arguments [6]. An 
advantage of eq. (|10p is that its coefficients are written in terms of observable quantities. Indeed, 
PQ + PQ is proportional to 1 + w, where w = PqI Pq is the equation of state of quintessence, which 
we will take here and in the following to be constant. The parameter M 4 is related to the speed of 
sound of quintessence, given by 

2= PQ+PQ nn 
S PQ+PQ+4M* ' 1 ' 



2 We assume that 4> is a monotonous function of t. 

3 The time variation of M 4 is expected to be even slower than Hubble, i.e. of order (1 + w)H, which is the 
typical time variation of pq. 
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As can be seen from this equation, absence of ghost - i.e., positiveness of the time kinetic-term in 
eq. (fTUj) - implies that c 2 s has the same sign as 1 + w [35l EE S]- In particular, for w < —1 one 
has c 2 . < 0, which signals the presence of gradient instabilities. As shown in d] stability can be 
guaranteed by the presence of higher derivative operators but requires that the speed of sound is 
extremely small, practically zero [lj. 

Regardless of the motivations expressed above on the stability of single field quintessence for 
w < — 1, in the following we will be interested in considering the limit c 2 — > 0, which is obtained 
when \Pq+Pq\ <C M 4 . We will see that what turns out to be physically relevant are the density and 
pressure perturbations on surfaces of constant (ft, i.e. of constant tt. These are the perturbations in 
the so-called velocity orthogonal gauge, and using eq. ((8|) they are given by 

6 Pq' ) = 6 PQ ~ PQ^ = {.PQ +PQ + 4M 4 )t^ , 5p { g ■°- ) = 5p Q - Pqtt = (p Q + p Q )n . (12) 
Indeed, c 2 defined in eq. (fTUj) can be written as [I] 

<* = d • < 13 > 



Thus, the pressure perturbation is suppressed with respect to the energy density perturbation by 
the smallness of the speed of sound. As we will see, in the limit c s — > this implies that pressure 
forces are negligible and quintessence follows geodesies, remaining comoving with the dark matter. 
In the limit c s — > 0, the energy density perturbation on velocity orthogonal slicing becomes 

5 P { q°- ] = AM 4 * . (14) 

Note that since -k ~ Htt, the difference between Spq'°'^ and 6pQ is negigible for small speed of sound, 
Spq' '^ ~ SpQ. All these conclusions hold independently of the value of Spq/Pq, provided that the 
effective theory described by action (|10p remains valid, i.e. for tt <C 1 [36, 4J. In particular, they hold 
also when perturbations in the energy density of quintessence become non-linear, i.e., for 5pQ 3> Pq- 

Gravitational perturbations can be straightforwardly included in the action (|10p as in [3]. As a 
warm-up exercise we will here, instead, study the evolution of quintessence in the spherical collapse 
solution. According to the spherical collapse model, the overdensity can be described as a closed 
FLRW universe with a scale factor R which is different from the one of the background a. This 
remains true also when we take into account quintessence with negligible speed of sound. Indeed, 
eq. (I13p shows that there is no pressure difference between the inside and the outside of the overdense 
region. Therefore, inside the overdensity we can describe quintessence using eq. Q, where the time 
evolution of the metric is described by the scale factor R and we thus replace a 3 by R? . 

With this new action, the second term proportional to tt can be integrated by parts and the 
coefficients of the linear part of the action rewritten in terms of pQ and pQ using eq. ©. However, 
now the linear part of the action does not cancel but can be written, using the background equation 
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of motion, as —5H(pQ + Pq)tt, where we have defined 5H = H m — H, with H m = R/R. We are thus 
left with a linear term in the action, due to the difference between the rates of expansion inside and 
outside the overdensity. After manipulations of the last two terms in eq. ([9]) , similarly to what was 
done to derive eq. (fTU|) . the action inside the overdensity becomei""' 



/ 



S= (fx R J 



\ {pQ+pQ + 4M 4 ) vt 2 - ^(p Q + p Q )^- + ~#(p Q +p q )tt 2 
3 

- 35H(p Q + Pq)tt - -5H(p Q + p Q ) it 2 



(15) 



Using that 7r ~ Htt, neglecting time variations of M and discarding terms suppressed when 
\pQ +PqI ^ (i- e -; i n the limit c 2 — > 0) the equation of motion of ir derived from this action reads 

R 35 H . _ _ . 

n + 3 R 7r = ~4M*^ Q+pQ >' ' 

As expected, the quintessence perturbation induced by the overdensity is proportional to 1 + w, i.e. it 
vanishes in the limit of the cosmological constant. Note that the source term on the right-hand side 
of this equation can be written as —3c 2 5H and is suppressed by the smallness of c 2 . This implies 
that, even for large over densities, i.e. 5H > H, variations of it due to the gravitational potential 
well are extremely small, 7r ~ c 2 , inside the regime of validity of the effective theory. Furthermore, 
this also implies that the difference in ir between the homogeneous and closed FLRW solutions is 
also tiny, Att ~ (? s H~ l . Thus, the quintessential scalar field practically lies on the same point of its 
potential. 

Equation (TTOj) can be written, using eq. (fT¥|) (and 5pQ ~ dpg' '^), as 

5p Q + 3^S PQ = -35H(p Q + p Q ) . (17) 

Note that, as 5H is always negative, the sign of 5pQ is the same as that of 1 + w. Remarkably, 
this implies that for w < — 1 dark matter halos accrete negative energy from quintessence, as 
was noticed at linear level in [9]. Combining eq. (|17p with the background continuity equation, 
Pq + 3H(pQ +Pq) = 0, we obtain 

PQ + ^{pq+Pq) = 0. (18) 

This equation describes the evolution of the energy density of quintessence with c 2 = inside a 
spherical overdensity dominated by dark matter. Note that the pressure perturbation is absent, as 
it is suppressed by <? s — > 0. Indeed, this equation differs from the description currently given in the 
literature for clustering dark energy. In particular, the analogue of this equation given in \27\ [29] 



4 We will not include in the action the contribution to y/—g coming from the curvature of the closed FLRW 
universe. Indeed, as it is time independent, it does not affect our discussion. 
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includes the pressure perturbation 5pQ = w5pQ. Including the pressure perturbation 5pQ leads to 
an incorrect description even in the linear regime, in contrast with eq. (|18p which does match the 
linear theory for small over densities. 

In the following two sections we will make this analysis more complete and derive all the equations 
necessary to describe the spherical collapse with quintessence. 



3 Spherical collapse in local coordinates 

As we did in the former section, the spherical collapse is usually treated using FLRW coordinates, 
as in the simplest cases the overdensity evolves as an independent closed universe. This somewhat 
obscures a crucial simplification of the problem, i.e. that the collapse of dark matter haloes occurs on 
scales much smaller than the Hubble radius. In this limit one can treat gravity as a small perturbation 
of Minkowski spaced As the dynamics of c s = quintessence is not completely intuitive, we want to 
make use of a coordinate system where this simplification is explicit; this will also make the dynamics 
of the other cases of spherical collapse clearer. We thus choose a coordinate system around a given 
point, such that the deviation of the metric from Minkowski is suppressed by H 2 r 2 , where r is 
the distance from the point, for any timely Notice we do not want to limit the validity of our 
approximation to times shorter than H -1 because this is also the typical time-scale of the evolution 
of a dark matter halo. These requirements define the so called Fermi coordinates. Note also that we 
are not taking any Newtonian limit: as we are interested in quintessence we cannot neglect pressure 
as source of gravity. 

A particular choice of Fermi coordinates are the so-called Fermi normal coordinates |38| where 
the deviation of the metric from Minkowski can be written as a Taylor expansion around the origin 
whose leading coefficients are components of the Riemann tensor. These are (with the convention 
of [39]) 

500 = -1- R iom\dX l x m + ... , (19) 
2 

1 



goi = 0- ^R i im \ 6 x l x m + ... , (20) 



9ij — &ij ^Riljrn\Q XX + . . . . (21) 

Here we will be interested only in spherically symmetric solutions. As Rou m vanishes because of 
rotational symmetry, goi must be of order higher than r 2 . Thus we can neglect it in the following 
discussion. Furthermore, rotational symmetry implies that the corrections to goo will be proportional 



5 For a recent use of this approximation in cosmology see [37]. 

6 In the spherically symmetric case, the range of validity of this approximation goes to zero close to the 
collapse singularity. However, this is not relevant because the singularity is anyway an artifact of the spherical 
symmetry. In the real case the curvature of space remains small and the halo reaches virialization. 



S 



to r while those to gtj will be proportional either to r or to x\Xj. It is possible to make a redefinition 
of the radial coordinate such as to get rid of the term X{Xj without affecting goo an d goi at 0(r 2 ). 
Note that in such a way we are using Fermi coordinates which are not of the normal form. In this 
case the metric can be written in the Newtonian gauge (not to be confused with the cosmological 
perturbation theory Newtonian gauge) form as 

ds 2 = -(l + 2$)dt 2 + (l-2^)dx 2 , (22) 

where and \& are proportional to r 2 . In this gauge the 00 component of the Einstein equation 
gives 

V 2 ^ = AirGp . (23) 
The part of the ij Einstein equation proportional to the identity gives 

6^ + 2V 2 ($-^) = 24vrGp. (24) 

As the typical time scale is of order Hubble, the \E' term is suppressed with respect to V 2 ^ by 
0(H 2 r 2 ) and can therefore be neglected. Thus, using eq. (f23|) we obtain 

V 2 $ = 4vrG(/) + 3p) . (25) 

As a first step, let us show how one can use these coordinates to describe an unperturbed FLRW 
solution with non-vanishing curvature. In isotropic comoving coordinates this metric is written as 



^-^♦•M'prgssi- (26) 



where K is the curvature parameter. With the change of coordinates r = t — \Hr 2 and y - 
-(1 + jH 2 r 2 ) [37], with a and H evaluated at t rather than at r, one gets at first order in H 2 r 2 , 

1 



ds' 



1-(H + H 2 )r 2 



dt z + 



l-^(H 2 + K/a 2 )r 2 



dx z , (27) 



which is indeed of the Fermi form (|22p , where the corrections from flat spacetime are given by 

$ = --(H + H 2 )r 2 , * = ^(H 2 + K/a 2 )r 2 . (28) 

Let us now show that the metric (|27p is a solution of the Einstein equations f)23[) and (|25p . In 
the coordinates (t,x l ), p[r) and p(r) are not space independent; however, their space dependence is 
suppressed by H 2 r 2 so that it can be neglected. With spherical symmetry, assuming regularity at 
the origin the two Einstein equations (|23p and (|25p are then solved by 

8-irGpr 2 . . 

* = — j , (29) 

and 

AttC r 2 

* = ^f(p + 3p)-. (30) 
9 



Comparing these expressions with ([25]) we recover the two Friedmann equations, respectively, 

2 K 8ttG 

H + — = -r-p , (31) 

and 



a AttG , , . 

- = 5~(p + 3p). 32 

a 6 

Note also that the traceless part of the ij Einstein equation, (d{dj — |<$yV 2 )(<& — x &) = 0, is trivially 
satisfied by the expressions above. Matter stays at fixed y in the original FLRW coordinates; 
therefore it moves in the new coordinates as x oc a, i.e. with velocity v = Hx. Finally, using this 
equality one can check that also the Oi component of the Einstein equation is satisfied. 

Let us now look at the dynamical equations for the fluid. The time component of the conservation 
of the energy- momentum tensor gives (see for example [40] ) the continuity equation, 

p + V- [(p + p)v\ =0 . (33) 

This equation is the same as in Minkowski spacetime as the gravitational corrections only induce 
terms suppressed by 0(H 2 r 2 ). When the velocity v is simply given by an unperturbed Hubble flow 
we obtain the standard conservation equation in expanding space, p + 2>H{p +p) = 0. 

The spatial component of the conservation of the energy-momentum tensor gives the Euler 
equation, 

V$ , (34) 



v + (v ■ V)v 



_ dp 
Vp + v- 



{p + p) 

where we have assumed that v <C c. At leading order in 0(H 2 r 2 ) gravitational perturbations 
enter only through the last term on the right-hand side of this equation. In the particular case 
of an isotropic and homogenous solution the first term on the right-hand side exactly cancels: as 
Vp = —pHx, the gradient of the pressure cancels with the term coming from its time dependence. 
This is not surprising as what matters in the Euler equation is the 4-dimensional gradient of pressure 
perpendicular to the fluid 4-velocity. In this case eq. ([34"]) reduces to 

?7+ (v ■ V)u = -V$ . (35) 

This equation is verified by the Hubble flow v = Hx since we get 

-x = -V$ , (36) 
a 



which is clearly satisfied by the explicit expression for $, eq. ([28]). 

We can now use these local coordinates to describe the spherical collapse in various cases, starting 
from the simplest. 
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• Dark matter only 

Let us take a spherically symmetric distribution of matter around the origin. As both the gravi- 
tational potentials $ and \& satisfy the Poisson equation, we do not need to know how the mass 
is radially distributed to solve for the gravitational background outside a given radius r. We just 
need the total mass inside the radius r. In particular (see figure [I]) if inside a given radius r out , a 
distribution contains as much matter as the unperturbed cosmological solution, from the outside 
it will look exactly as the unperturbed background. This implies that we can smoothly glue this 
solution at r = r ou t to the cosmological background, and that the latter will not be affected by the 
gravitational collapse inside. This is of course a linearized version of Birkhoff 's theorem in General 
Relativity. 




Figure 1: Spherical collapse 

Conversely, the solution inside a given radius is not affected by what happens outside. In 
particular, if we assume a homogeneous initial condition inside a radius r; n (with r; n < r out ), this 
central region will evolve as if these homogeneous initial conditions were extended outside, i.e. as 
a complete FLRW solution |17j . The central overdense region will remain exactly homogeneous, 
reaching maximum expansion and then collapsing. 

Without further assumptions, the evolution of the layer r; n < r < r out does not enjoy particular 
simplifications and its evolution must be computed as a function of the initial profile. In any case, 
this is usually irrelevant as we are only interested in the fate of the r < r; n region. If we assume 
that the layer r ia < r < r out is empty, the Poisson equations ([23]) and ([25]) give 1/r solutions for the 
potentials, like for a source localized at the origin. This is the linearization of the exact Schwarzschild 
solution. 
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• Dark matter and a cosmological constant 

The considerations above also hold when we include a cosmological constant. Although $ and ^ 
now solve different equations (because p\ 7^ 0), they are both of the Poisson form. Thus, we still 
have an unperturbed evolution outside r out if the total matter inside matches the background value. 
Assuming initial homogeneity, the central region r < rj n will evolve like a complete FLRW universe 
|23j . Although now pressure does not vanish, it just comes from the cosmological constant which 
does not define a preferred frame and is therefore comoving with dark matter both inside and outside 
the overdensity. 

• Non-clustering quintessence: c s = 1 

When quintessence has a speed of sound c s = 1, it does not effectively cluster but it keeps on 
following the cosmological background solution, irrespective of the dark matter clustering [24] . 

As before, outside r out there is an unperturbed cosmological background. What is new now is 
that the central region r < r- m does not behave as a complete FLRW solution, even if we start with 
a homogeneous overdensity. Indeed, quintessence and dark matter do not have a common velocity: 
while dark matter slows down and eventually starts collapsing, quintessence keeps following the 
external Hubble flow vq = H ou tx. Note that, on the other hand, in the cosmological constant case 
one cannot define a dark energy 4-velocity as its energy-momentum tensor is proportional to the 
metric. To study the evolution of the dark matter overdensity one must use the Euler equation (|36p . 
Here, what defines the velocity flow of dark matter is the effective "scale factor" R, v m = R/R x. 
This yields 

£ = _V$. (37) 

R 

Using the explicit solution (f3"Uj) for $, this equation becomes [M] 

R 4ttG . _ , . . 

^ = ^-{Pm + PQ + 3PQ) , (38) 

where we have separated the contribution of dark matter and quintessence to Notice that, 
although this equation looks like one of the Friedmann equations, the dynamics of R is not the same 
as for a FLRW universe. Indeed, p m evolves following the scale factor R, while the quintessence 
follows the external scale factor a. In a FLRW universe, from eq. f|38[) together with the continuity 
equation one can derive the first Friedmann equation, (R/R) 2 = ^j^-p — K/R 2 . Here, as the different 
components follow different scale factors, this is not longer possible and the first Friedmann equation 
does not hold. 

• Clustering quintessence: c s = 

Let us now move to the subject of this paper. We want to show that in the limit of vanishing speed of 
sound quintessence remains comoving everywhere with dark matter. In particular, this implies that 
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in the region r < r- m quintessence follows dark matter in the collapse and the overdensity behaves 
as an exact FLRW solution so that, contrary to the c s = 1 case, also the first Friedmann equation 
holds. The fact that quintessence remains comoving with dark matter can be understood both by 
using the fluid equations or directly from the scalar field equation of motion. 

In the fluid language the dynamics is described by the Euler equation (|34p . In general, in the 
presence of sizable pressure gradients a fluid does not remain comoving with dark matter, i.e. it 
does not follow geodesies. Since quintessence has a sizable pressure, the fact that it moves following 
geodesies may be unexpected but it is obtained in the limit c s — > 0. This can be easily seen by 
rewriting the Euler equation for quintessence in covariant form as 

U » V v = - 1 U"> + u»u°)V aPQ , (39) 
[PQ+PQ) 

where vP" is the quintessence 4-velocity. When the right-hand side of this equation vanishes, the 4- 
velocity solves the geodesic equation. Notice that the pressure gradient is multiplied by the projector 
perpendicular to the fluid 4-velocity. This is the same as projecting on surfaces of constant 4> an d it 
is equivalent to a gradient of the velocity-orthogonal pressure perturbation that appears in equation 
(|12p . which involves only tt, and not tt. By eq. (|13p this is negligible in the limit c s — > and thus 
the right-hand side of (|39p vanishes. 

This result is even clearer in the scalar field language. Taking the derivative of the equation 
defining the quantity X in ([T]) , 

d y {d^(j)d^) = -d u X , (40) 
and writing it in terms of the 4-velocity v} 1 = —d^(pj\fX we have 

2?/ t VxV M (\/Xu ! = -d u X , (41) 

and therefore 



t^VX = -i^W + u v u a )d a X . (42) 



1 

2X 

Equation ([39]) is recovered using eq. ([5]) and taking into account that d a P((j),X) = dP/dcp ■ d a (j) + 
dP/dX ■ d a X and that the first term vanishes when multiplied by the projector orthogonal to 
vP" . From this we clearly see that what matters is only the gradient of the pressure on <p = const 
hypersufaces. This vanishes in the limit c s — > and thus we have geodesic motion. We stress that, 
although quintessence with c s = follows geodesies, its dynamics is quite different from dark matter. 
Pressure does not accelerate the quintessence 4-velocity but it does affect the energy conservation 
equation (|33|) . Moreover, quintessence does not enjoy a conserved current, while dark matter particle 
number is conserved; this is related to the absence of the shift symmetry (p —> <p + c in the scalar 
field Lagrangian (see for example |41|). 

As discussed in section [21 the different dark matter evolution inside and outside the overdensity 
changes the quintessence solution by a very tiny amount A-zr ~ (? s H~ l : the quintessence field sits 
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at the same position along its potential, (ft = (ft(t), apart from negligible c 2 corrections. Notice 
that this was derived using two different Friedmann coordinate systems, one following dark matter 
inside the overdensity and one following the unperturbed Hubble flow outside. Thus, in reality we 
have two solutions (ft = 4>(ti n ) and (ft = 0(t O ut) respectively. Once these two solutions are written 
in the same local coordinates (|27p . the solution for (ft becomes (ft = (ft(t — ^Hi n r 2 ) for r < r m and 
(ft = (ft{t — ^H out r 2 ) for r > r out . This implies that in these coordinates it has to jump in the layer 
between the two regions, by an amount Att ~ ^r 2 5H, and that this jump is not suppressed by 
c 2 . One may expect that the scalar field would "react" to this gradient between the inside and 
the outside. However, this does not happen in the limit c s — > as the spatial kinetic term is very 
suppressed. Let us see this explicitly. 

To study the scalar field equations in the local coordinates, one can start by writing the equations 
in Minkowski space and then check a posteriori that the deviation of the metric from flat space only 
gives relative corrections 0(H 2 r 2 ). The evolution equation for (ft, eq. ([2]), reads in Minkowski space 

- d t (P x <j>) + diiPxdicft) = -\p^ . (43) 

If we try a comoving solution of the form (ft = (ft(t — ^Hr 2 ) we end up with the standard FLRW 
equation for (ft, inclusive of the friction term, 

- d t (P,x$) - 3HP X $ = -~P |0 . (44) 

Metric fluctuations give only a correction to this equation of order 0{H 2 r 2 ). As we discussed, the 
two homogeneous solutions for r < r- m and r > r ou t are different so that we expect gradient terms 
to smooth out the initial top-hat profile. To estimate the thickness L of the layer over which the 
smoothing takes place, we can study perturbations around a top-hat profile and require the spatial 
and time kinetic term of the perturbation it in eq. (jlOp to be comparable, 

MV~(pq+Pq)^. (45) 

Using that tt ~ Hit, from this comparison we obtain L ~ |c s |-ff -1 . This makes perfect sense: our 
top-hat profiles are smoothed out over a distance comparable to the sound horizonQ 

In conclusion, the solutions outside and inside the overdensity are exact FLRW with quintessence 
comoving with dark matter. Gradient terms will smooth out this solutions on scales of order of the 
sound horizon, which vanishes for c s — > 0. This discussion also tells us that taking c s = will be 
correct only for objects which are much bigger than the sound horizon |c s |-?/ _1 . In the opposite 
limit of an object which is much smaller than the sound horizon, one can treat quintessence as 
unperturbed as discussed above in the c 2 = 1 case. For example, if one is interested in objects larger 
than 1 Mpc, one can neglect the speed of sound as long as \c s \ < 10~ 4 . 

7 It is straightforward to check that this estimate is not altered by the higher derivative operators that are 
required for stability when w < — 1 [5]- 
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4 Solving the spherical collapse 



In this section we derive the equations for the spherical collapse of dark matter in the presence of 
quintessence with vanishing speed of sound and we compute their solutions numerically. 

• The background universe 

The background is described by a flat FLRW metric with scale factor satisfying the Friedmann 
equation, 

a J 3 

where p m and pQ are the background energy density of dark matter and quintessence, respectively. 
For later purposes, we express p m and pQ in terms of the fractional abundance of dark matter J7 m , 

- — 3 ^ 2 <~> - 1 - _ 

Pm = -Mm , PQ = —T, Pm ■ (47) 

Dark matter redshifts with the expansion as the physical volume, p m oc a -3 , while the energy density 
of quintessence scales as pq oc a ' W K The dark matter contribution to the critical density f2 m 
can be written as a function of its value today, O mi o, and x, the scale factor normalized to unity 
today (at t = to), 

x = a/a Q ■ (48) 



This yields 



HJr) + . (19) 



-1 



Equation (|47p can be then rewritten as 

87TG X 6 "m,0 
where the second equation follows from (|49p . Furthermore, rescaling the time variable by defining 

V = \AW , (51) 

one can rewrite the Friedmann equation as 

^ = (xfl^x))- 1 / 2 . (52) 

The initial condition for x can be imposed at some small initial time rji during matter dominance, 
Xi = (3r/j/2) 2//3 . Then, eqs. (|49p and ()52p completely describe the background evolution of the metric 
and energy-momentum tensors. 
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• The linear evolution 

Before studying the collapsing spherical overdensity we derive the evolution equations of perturba- 
tions of dark matter and quintessence in the linear regime. As we consider scales much smaller than 
the Hubble radius, the gauge dependence is not important. We will thus perturb the continuity and 
Euler equations in local coordinates, eqs. (|3"3"j) and adding small inhomogeneous perturbations 
6(t, x) and u(t, x) to the homogeneous energy density and Hubble flow velocity, 

p = p(l + d), v = Hx + u. (53) 

Let us start from the dark matter. Perturbing at linear order eq. (|33p with p m = yields 

d -+\ 

— + Hx -VU m = -V • u (local coords) , (54) 

where we have specified that we are describing perturbations using local spatial coordinates x. On 
the other hand, on the left-hand side of this equation one recognizes the time derivative at fixed 
comoving coordinates y = x/a(t), i.e., 

(55) 

Indeed, here we are interested in describing the evolution of an overdensity of dark matter contained 
in a comoving volume. Thus, we describe 5 m and u as a function of the comoving coordinates, which 
simply gives 

5 m + -V • u = . (56) 
a 

To close this equation we need the evolution of the dark matter peculiar velocity u. This can 
be obtained by perturbing at linear order the Euler equation (|35p . Using comoving coordinates the 
perturbed Euler equation becomes 

il+Hu+ -W$ = , (57) 
a 

where <5<J> is the perturbation of the Newtonian potential, 

8$ = $ + ]-(H + H 2 )r 2 . (58) 



Equations (|56p - (|58p have been derived for instance in [52] in the context of Newtonian mechanics 
described with expanding coordinates, for a pressureless fluid in the presence of vacuum energy. 
Here the Poisson equation for 5& is sourced by both dark matter and quintessence perturbations, 

-lv 2 5$ = 4irG(p m 5 m + p q 5q) , (59) 

where we have used that 5pQ = 0. The final step is to eliminate the peculiar velocity by subtracting 
the divergence of eq. (|57p from the time derivative of eq. (|56p. With the Poisson equation (|59p we 
obtain 

'5 m + 2H5 m = 4TrG(p m 5 m + pq5q) . (60) 



16 



For quintessence we perturb the continuity equation (f3"3"|) which gives, in comoving coordinates, 
using Sp Q = 0, 

6q-3Hw5q + (1 + w)-V -u = . (61) 



To eliminate the divergence of the peculiar velocity we can use eq. (j56[) taking quintessence to be 
comoving with dark matter. Indeed, as explained above, both the dark matter and quintessence 
follow geodesies and are dragged by the same potential well and the growing mode of their velocities 
is the same. Thus 

5 Q - 3Hw5 Q = (1 + w)8 m . (62) 

In matter dominance, when 5 m oc a, the solution of this equation is [4J 

l + w „ 



(63) 



Note that the denominator on the right-hand side further suppresses the perturbation of quintessence 
with respect to the naive 1 + w estimate. 

In terms of the dimensionless variables x and rj, respectively defined in eqs. (|48p and (|5ip . 
equations ([60]) and (j62|) rewritten as 

d 2 5 m 2dxd5 m _ 3 / 1 - fl mfi _ 3w 



where we have used eq. (|50j) . and 

^9. - 1**E. W 8 _ | d£ m 
dr/ a; d?7 d?? 

The initial conditions are set in terms of the initial dark matter density contrast 8 m ,i- I* 1 matter 
dominance 5 m = H5 m i.e., dS m /dr)\i = 25 m) j/(3ryj), while the value of Sq^ is fixed by through 
equation ([63]) . 



• The spherical overdensity 

We now study the evolution of a spherical homogeneous overdensity of radius R in a FLRW back- 
ground that satisfies the Friedmann equation (|46p . We denote the energy densities of dark matter 
and quintessence inside the collapsing ball by p m and pq, respectively. Since dark matter is pres- 
sureless p m = and since quintessence pressure perturbation is negligible, 5pQ <C §Pq, we can take 
quintessence pressure to be the unperturbed one Pq. 

In local coordinates, the evolution of the scale factor R is described by the Euler equation (|36p . 
Using the appropriate scale factor - i.e., R instead of a - and replacing the potential $ using eq. pOj) . 
the divergence of this equation can be written as 

R 4vrG . , 

P = ^~ {Pm + PQ + 3p Q ) . (66) 
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(Note that for a non-clustering quintessence the equation for R is the same with pq replaced by 
PQ M ) 

For the evolution equations for p m and pq we use the continuity equation ([33]) . Inside the ball 
this reads, for dark matter, 

Pm + 3 ^Pm = , (67) 



whose solution is simply 



For dark energy eq. ([33]) becomes 



' R 3 ' 



PQ+^(PQ+PQ) = , (69) 
which can be rewritten in terms of the nonlinear density contrast Aq = Pq/Pq — 1 as 

A Q + 3|a q - 3^(1 + w)A Q + 3(1 + w) ^| - = . (70) 

To solve eqs. (|66p and ()70p numerically it is convenient to use y, the radius of the ball normalized 
to unity at the initial time, 

y = R/Ri , (71) 

and change a and t to the dimensionless variables x, r]. Using eq. (|50p . eq. (]66p can be rewritten as 



l 

dr/ 2 2 



1 + m j 1 . 1 f^m,o V 

+ (1 + 3tu + A c 



0, (72) 



where we have used eq. (|68|) and that in the linear regime, where the initial conditions are set, 
Pm/Pm\i = 1 + 5 m>i . Equation (JZOJ) yields 

dAn . N /ldy l + AndxX A dw , . 

y +3(l + u>) --^ — 3~ +3— r = - ( 73 ) 



dr] \y dr] x dr] J y dr] 

As initial conditions we have yi = 1 by definition; the expansion rate of a collapsing sphere with 
dark matter only and in the linear regime can be written as 



1=1 1 -I s -)' < 74 > 

which fixes the first derivative of y, dy/dr]\i = 2(1 — 5 mj j/3)/(3r/j). For the dark energy perturbation 
we use that Aq^ is linear at early times, Aq^ = Sqj, and thus is fixed in terms of S mi i by eq. ([63]) . 

By solving numerically eq. (|52|) for the background evolution described by x and plugging the 
result into the coupled eqs. ([72]) and ([73]) . one can compute the evolution of R as a function of time 
t, from the initial time ti to the time of collapse t c . The evolution of R is shown in figure [2] for 
four different models: CDM only, ACDM, c s = 1 quintessence and c s = quintessence in the cases 
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Figure 2: Thick lines: time evolution of the radius for a spherical collapse. Thin lines: time 
evolution following the linearized solutions. The quintessence equation of state is w = —0.7 (above) 
and w = —1.3 (below). Starting with the same overdensity, a model with CDM only is the first to 
collapse. In the upper figure ACDM collapses before the quintessence models as dark energy with 
w = —0.7 is more important in the past. The situation is reversed for w = —1.3. For w = —0.7 
the c s = quintessence collapses before c s = 1 as positive energy clusters together with dark matter. 
For w = —1.3 the situation is reversed as negative energy clusters and hinders the collapse. Note 
that quintessence models with c\ = 1 and w < — 1 are plagued by ghost instabilities and are thus 
very pathological on short scales. In this figure and in the following ones we study this case only for 
comparison with the c\ = case. 

w = —0.7 and w = —1.3. We stress that quintessence models with c 2 s = 1 and w < — 1 are plagued 
by ghost instabilities and are thus very pathological on short scales. We study them here only for 
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comparison with the c 2 = case. We asume £l m ,o and Oq,o (or ^a) to be the WMAP5 best fit values 
[44] . We have taken = 3 • 10 -4 as initial dark matter overdensity at rji = 1CP 6 . As expected, 
since the cosmological constant and the quintessence slow down the evolution of R, the collapse is 
faster in the pure CDM model. This effect takes place earlier for w > —1 as in this case quintessence 
is more important in the past than the cosmological constant. Thus, for w = —0.7 the collapse 
happens later. On the contrary, for a quintessence with w = —1.3 collapse takes place earlier. For 
w = —0.7 the collapse is enhanced by the quintessence perturbations and it takes place faster when 
c 2 = 0. The opposite happens for w < —1, as in this case negative energy clusters, hindering the 
collapse (see eq. (fT7|l ). 

In general, the time of collapse depends on the value of the initial dark matter overdensity. 
This is shown in figure El where the redshift of collapse z c is plotted as a function of the initial 
density constrast 5 m ^ at the same initial time. As expected, larger ovedensities collapse earlier, at 
higher redshift. For large enough overdensities - and early enough collapse - the redshift of collapse 
becomes the same for all four different models, because the cosmological constant or the quintessence 
remain subdominant during the whole process. As expected, quintessence with w > — 1 requires a 
larger initial overdensity to collapse and in this case quintessence perturbations (c 2 = 0) help the 
collapse. The opposite happens for w < — 1. 

An important quantity to compute in order to derive the mass function is the critical density 
contrast 5 C , i.e., the density contrast in the linear theory computed at the time when the spherical 
collapse solution reaches the singularity. Thus, we numerically solve the linear evolution equations 
for 5 m and 5q, eqs. ([64^) and ([65]) . and we take S c to be 5 m at the time of collapse. In the standard 
CDM scenario 5 C is given by the well-known number 1.686 [17] independently of the redshift of 
collapse z c . However, in the presence of a cosmological constant or quintessence, 5 C depends on 
the redshift of collapse. Indeed, as the relative abundance of dark matter and dark energy changes 
with time, the dynamics of the spherical collapse depends on when it takes place. This is shown in 
figure [J where we plot 5 C as a function of z c . This result generalizes to quintessence with c 2 = 
the standard results obtained for CDM [17] . ACDM [23] and smooth quintessence [2j. As expected, 
if the collapse takes place early, when the cosmological constant or quintessence are not important, 
the critical density 5 C will be the same as for CDM. The cosmological constant decreases the value of 
5 C and quintessence with w > —1, becoming important earlier, decreases it even more. Quintessence 
perturbations enhance S c if w > — 1. This effect is very mild for w < — 1 because quintessence 
becomes important only at very late time. 

The change of threshold is very small (< 0.5%) [28] in all the cases and it is easy to understand 
why. Let us compare for example a universe with CDM only with a ACDM one and let us focus on 
objects that collapse at a given redshift, say z = 0. The initial overdensity must be rather bigger 
in the ACDM case to overcome the acceleration induced by A. But the threshold 5 C is obtained 
evolving this initial value with the linear transfer functions and this will suppress the ACDM value, 
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Figure 3: Redshift of collapse as a function of the initial overdensity. In the upper figure quintessence 
models have w = —0.7, while they have w = —1.3 in the lower one. The behavior follows that 
explained in figure 

exactly for the same reason which required it to be bigger in the first place. In other words the only 
effect comes from the difference between the linear and non-linear evolution and this causes only a 
small suppression with respect to the CDM case. 

5 The mass function of dark matter halos 

We are now ready to discuss the predictions for the mass function using the Press-Schechter formalism 
|18} 119] , We will first concentrate on the mass of dark matter, leaving aside for the moment the 
contribution of the quintessence mass to the halos. The volume density of dark matter halos of mass 
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Figure 4: Linear overdensity at collapse as a function of the redshift of collapse. In the upper figure 
the quintessence models have w = —0.7, while they have w = —1.3 in the lower one. 



M is given by 



^PS/jit \ [5 p m 5 c (z) d\oga M 



6: 



2/ 



(75) 



dM v ' ' V TrM 2 D(z)a AI dlogM 
Here cr^ is the smoothed variance of the density field which we define with a sharp cut-off in real 



space 



1 r 00 Att 
a 2 M = ^j^ dkk 2 \W(kR)\ 2 P m (k) with M = —R i p m , (76) 

where P m is the matter power spectrum and W(kR) = 3(sin kR — kRcos kR)/(kR) 3 is the Fourier 
representation of the top-hat window function in real space. Note that the Press-Schechter mass 
function (|75p can be rigorously derived only using a sharp filter in Fourier space. Thus its use 
with a sharp filter in real space is just an approximation; for corrections to this approximation see 
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for example [2T]. Notice that the redshift dependence of the threshold 5 c (z) only comes from the 
spherical collapse dynamics discussed in the previous Section and does not include the growth of the 
matter power spectrum, which is separately taken into account by the linear growth function D(z). 

The linear matter power spectrum is very similar in the cases c s = 1 and c s = 0. The difference 
comes from the contribution to the Poisson equation of quintessence perturbations as shown in 
eq. (|64|) . The effect is independent of k and it is thus equivalent to a change in the growth function 
D{z). The change in the growth function can be easily calculated as a function of z by solving 
eqs (|64p and (|65p . The result is shown in fig. [5j Given that quintessence becomes relevant only 
recently and that perturbations are suppressed by 1 + w the effect on dark matter does not exceed 
the percent level; this result is somewhat smaller than a naive estimate one can get by comparing 
the two contributions on the right-hand side of eq. (|64|) . As one can see in fig. for 1 + w > setting 
the speed of sound of quintessence to zero fosters the clustering and the dark matter spectrum is 
slightly enhanced; for 1 + w < quintessence clustering has negative energy and the dark matter 
spectrum is suppressed. The size of the effect is smaller for 1 + w < as quintessence becomes 
relevant only very recently. 
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Figure 5: Ratio between the growth functions D(z) for c s = and c s = 1 as a function of the redshift 
z. 

For the calculation of the mass function, we do not need only the growth rate, but the complete 
matter power spectrum. We use for this the publicly available code CAMB |45j . which allows to set 
to zero the speed of sound of quintessence. Apart from w all the other cosmological parameters are 
set to the WMAP5 best fit values [2]. We have checked that the effect of setting c s = instead of 
c s = 1 in the code is compatible with what we got in fig. [SJ In figure we show the matter power 
spectrum with c s = and c s = 1 for two different values of w. The speed of sound gives a small 
effect, much smaller than the modification of the growth rate induced by the different background: 



23 



for w > — 1 quintessence becomes relevant before and suppresses the spectrum in comparison with 
ACDM. The opposite effect is obtained in the case w < — 1. 




Figure 6: Matter power spectrum. The two lines with w = —1.3 are almost superimposed. 

The equation of state w and the speed of sound of quintessence enter in the mass function 
eq. ([75]) modifying the growth function D{z) and the threshold for collapse 5 c (z). Notice that only 
the combination 5 c {z) / D(z) enters in the Press-Schechter formula. As we discussed in the previous 
Section, the change in the threshold is very suppressed and much smaller than the correction of the 
growth function. 

We want to focus on the effect of clustering quintessence on the collapse with respect to the 
case when quintessence remains unperturbed. This effect can be estimated taking the ratio of the 
Press-Schechter mass functions in the two cases, 

dn PS /dM(w,c s = 0) 

dn PS /dM(w,c s = 1) ' 1 ' 

plotted in figure UJ We see that the effect is quite small: the ratio becomes large at high mass as a 
consequence of the exponential dependence of the number density on the mass. We do not dwell on 
the measurability of this small effect because, as we will discuss in the next section, the additional 
contribution to the mass of the halo coming from the clustered quintessence will give a comparable 
change in the mass function. 

It is well known that the Press-Schechter formula does not fit in detail the mass function obtained 
by numerical simulations. A better fit, motivated by the ellipsoidal collapse model, is given by the 
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Figure 7: Ratio of the Press-Schechter mass function for c s = and c s = 1 at z = (above) and 
z = 1 (below). 



Sheth-Tormen mass function 1 20 
dn S T 



dM 



(M,z) 



12a , 
— A 

7T 



1 + 



ab 2 n 



D(z)*a 2 M 
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exp 



2D(z)W M 



(78) 



with a = 0.707, A = 0.322184 and p = 0.3. Since, as discussed, the dependence of the threshold 
5 C on the cosmology is very mild in all cases, 5 C is usually taken to be z independent and equal 
to the EdS value, 5 C = 1.686. Notice that in this way the function is "universal" in the sense 
that the dependence on the cosmological parameters and redshift is only through the smoothed 
linear density field D{z)gm- It is reasonable to expect that the Sheth-Tormen formula gives a 
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good description of the mass function in the case of non clustering quintessence. In this case the 
only effect of quintessence is through the time dependence of the background and its effect on the 
growth function; this is not qualitatively different from the case of ACDM. In other words we expect 
"universality" to hold also on this case. On the other hand, we can estimate the effect of clustering 
quintessence using the ratio described above: in taking the ratio we expect that the shortcomings 
of the Press-Schechter prescription will partially cancel. Therefore, in figure [5] we plot 

dn _ dn S T , dn PS /dM(w,c s = 0) 

dM ~ dM [W,Cs ~ ' ' dn PS /dM{w,c s = 1) • [ ' 

As expected the main effect is at low redshift and high masses. 



6 Quintessence contribution to the halo mass 

So far we have been interested in the contribution of dark matter to the halo mass function. Given 
that quintessence with vanishing speed of sound participates in the collapse, one may wonder whether 
quintessence will contribute a sizable amount of mass to the dark matter halo. After all, most of 
the measurements will be sensitive to the total mass of the object, and not only to the fraction of it 
associated with dark matter. 

The quintessence contribution to the halo mass can be defined as 

M Q = [ d 3 x 5 PQ , (80) 
Jm 

where the integral is extended over the whole dark matter overdensity. If we stick to the spherical 
collapse model this just reduces to (4tt/3)R?5pq. Of course it makes sense to interpret this expression 
as a contribution to the halo mass, only if it stays practically constant over the time scales of interest, 
i.e. a Hubble time. In the spherically symmetric case 5pQ follows equation ([Tj 



+ 3^Sp Q = -36H(p Q +p Q ). (81) 

When \5pq\ 3> |1 + w\p~Q, the right-hand side is negligible. In this limit 5pQ redshifts as matter 
so that the integral ()80p becomes constant as can be seen in figure (9) We expect this condition 
to be marginally satisfied at turn-around when 5pQ ~ (1 + w)pQ. This allows us to estimate the 
quintessence contribution to the halo mass, 

Mq , Mo 

+ (82) 

although for a precise estimate one cannot neglect the evolution of the quintessence mass after 
turn-around. Notice also that quintessence with 1 + w < contributes with a negative mass. 

Of course this is only the prediction of the idealized spherical collapse solution. In reality dark 
matter halos virialize with an overdensity ~ 200 times larger than the background. What happens 
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Figure 8: Mass function calculated using eq. ()79p for z = (above) and z = 1 (below). 

to quintessence while dark matter virializes? Quintessence is exactly comoving with dark matter; 
eventually dark matter reaches shell crossing and the velocity field ceases to be single- valued. This 
corresponds to the formation of cusps in the quintessence field, similarly to what discussed in the 
ghost condensate case [36] and more recently in the context of Hof ava-Lifshift gravity \4Q\ Wf\ . The 
dynamics of the cusps will depend on higher derivative operators and possibly on the UV completion 
In any case, the dynamics of quintessence in this phase is very complicated [36] 



of the theory 



3 For positive 1 + w one has c 2 s > 0. In this case, one would naively expect that quintessence remains 



27 



-0.05 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

'1 




Figure 9: Spherical collapse for w = —0.9 (above) and w = —1.1 (below). Solid blue line: the 
quintessence contribution to the halo mass (4tt/3)R 3 5pq, normalized to the constant dark matter 
mass, as a function of time. Note that it becomes constant at late time. Red dotted line: the same, 
but for the background energy density pq. Orange dashed: the same for (1 + w)pQ. 

and its treatment is beyond the scope of this paper. Let us assume however that the cusps are 
somehow regularized and try to draw some general conclusion that is independent of the details of 
virializationo 

As 5pQ ~ (1 + w)pQ at turn-around, when matter starts virialization the inequality \5pq\ 3> 
|1 + w\pQ is satisfied with rather good approximation: as the virial radius is approximately half of 
the turn-around radius, 5pQ has grown by approximately 8 times from the beginning of the collapse 
while pQ has remained approximately constant. The spherical collapse solution indicates that the 

smooth on very short scales thus preventing the formation of caustics. However, in our case the velocity of 

the quintessence fluid (which is the same as the dark matter velocity) exceeds the speed of sound, i.e. it is 

"supersonic" . In this case sound waves are too slow to prevent the formation of caustics. 

9 We are implicitly assuming that the process responsible for smoothing the cusps does not lead to an 

energy loss to infinity. For instance, this happens if the smoothing excites new relativistic degrees of freedom 

which are radiated away. 
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quintessence mass is with good approximation constant at virialization. In the real case we expect 
the variation of the quintessence mass to be even smaller; the conservation of the stress-energy tensor 
tells us that the mass variation is related to the energy flux across a surface around the object, 

M = jdSiT l Q . (83) 

As during virialization the velocity field of quintessence will cease to be radial, we expect this 
integral to be suppressed with respect to the spherically symmetric case. We conclude that, as 
the flux integral is negligible, the mass associated to quintessence stays constant independently of 
the details of virialization. Thus, a good estimate of the quintessence mass can be obtained from 
the spherical collapse model evaluating the quintessence contribution at the virialization radius, see 
figure [TUJ The ratio between the virialization and the turn-around radii is taken to be the same as 
in ACDM [23J@ 







\ 




> 


0.04 
0.02 


\ 

\ 

\ 

\ 

» \ 

\ N 


w=-0.8 ] 

w=-1.1 

w=-1.2 


> 








ii 

3 


0.00 
-0.02 


i 



1 2 3 4 5 



z 



Figure 10: The mass contribution of quintessence to the total halo mass, calculated from the spherical 
collapse solution when the radius reaches the virialization value. 

Let us see how we can take into account the additional quintessence mass in the prediction for 
the mass function. A rigorous treatment would be quite challenging: the formation of an object of 
mass M at redshift z should be accompanied by an extra mass M — > M[l + e(z)], where e(z) = 
^Q,vir/M m ^i r is the quintessence to dark matter mass ratio at virialization. Then one should follow 
this extra mass as the halo merges to form larger objects which in turn accrete extra quintessence 
as they form. However e(z) is important only at low redshifts, so that we expect the main effect to 
be on large objects that formed very recently, the ones on the exponential tail of the mass function. 

10 Quintessence clustering will modify the virial radius with respect to ACDM with corrections 0(\1 + w\). 
The effect of these corrections on the amount of clustered quintessence is 0(\1 + w\ 2 ) and can therefore be 
safely neglected. 
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Let us see how we can estimate the effect in this region. At any epoch, the largest objects are mostly 
created with a negligible rate of destruction through merging to form larger objects. Therefore the 
formation rate for these objects can be accurately approximated by 

-!**§-(M,z). (84) 
dzd\ogM K ' ; y ' 

As the effect of quintessence is to rescale the mass of the object as it forms, it is more convenient 
to use the mass function per logarithmic mass interval. Using this expression, when positive, as an 
approximation for the formation rate gives 

-{M, z) = — — e(z)), z) . (85) 



dz dlogM v ' ' dzdlogM 

Expanding the right-hand side of this expression and integrating it over the redshift enables us to 
take into account the extra mass associated with quintessence accreted from an inizial redshift Zi, 



dn PS , m +Q (M , dnps 
dlogM [M,Z) dlogM [M,Z) 

d dnps 



+ / dze(z) 



d d dnps 
' dlog Mdz dlogM 1 ,Z) 



dzdlogM 



(M,z] 



(86) 



where 6{x) is the Heaviside theta function. 

In figure [TT1 we plot the mass functions including the quintessence mass contribution, using the 
new Press-Schechter mass formula eq. (|86p into eq. (|79p . Notice that the effect of quintessence mass 
is to bring the c s = lines closer to the ACDM one. To better visualize this effect, in figure [12] 
we plot the ratio between the c s = and c s = 1 case. For w > —1 (w < —1), setting to zero the 
speed of sound of quintessence not only does it enhance (diminish) the clustering of dark matter as 
discussed in the previous Sections, but it also adds positive (negative) mass to the halo: the two 
effects therefore pile up and the second is quantitatively slightly dominant. The sum of the two 
effects is rather large: for values of w still compatible with the present data and for large masses the 
difference between the predictions of the c s = and the c s = 1 cases is of order one. 

We stress that the new mass function is not universal in the sense that there is an explicit 
red-shift dependence besides the one implicit in the growth of ctr This z dependence is quite 
remarkable: the fact that the modification takes place only at very low red-shift is quite distinctive 
and a clear link of the effect with the onset of acceleration. 



n The universality of the mass function, even in the presence of scale independent non-Gaussianity, has 
been tested with good accuracy in iV-body simulations [35]. 
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Figure 11: Mass function for z = (above) and z = 1 (below), including the quintessence mass 
contribution, calculated using eq. (|86p and (|79p . 

7 Three contributions to the mass 

In the previous section we saw that a distinctive signature of quintessence with c s = is the 
extra contribution of quintessence to the mass of virialized objects Although a detailed study 

12 In this section we will assume this mass to be positive, which is the case for 1 + w > 0. Similar consider- 
ations apply to the 1 + w < case. 
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Figure 12: Ratio of the Press- Schechter mass function, including the quintessence contribution 
eq. ([86]) . for c s = and c s = 1 at z = (above) and z = 1 (below). 

of how to distinguish this extra contribution goes beyond the scope of this paper, few remarks are 
in order. Let us focus on clusters: as discussed in the previous section, the mass of these large 
objects is significantly affected by the quintessence mass. Moreover, clusters are mostly dominated 
by gravitational physics (for an introduction to the subject see [12]). If we neglect quintessence for 
the moment, a cluster is characterized by its baryon mass, mostly in the form of gas, and by the dark 
matter mass. For sufficiently large clusters, in first approximation one expects the ratio between 
these two components to be close to the cosmological baryon to dark matter ratio |50j . 

The various techniques to study clusters have different sensitivities to the two components (for a 
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review on cluster detection tecniques see [51J). Weak lensing and galaxy velocity measurements are 
sensitive to the gravitational mass of the object, i.e. to the sum of the dark matter and baryon mass. 
On the other hand the optical richness, a measure of the light emitted by galaxies inside the cluster, 
is a probe of the baryon fraction of the mass. X-ray and Sunyaev-Zeldovich (SZ) measurements 
are sensitive both to the dark matter component and to the baryonic one. While these effects are 
caused by the gas in the cluster, the temperature and the equilibrium of the intergalactic medium 
depends on the dark matter distribution. Indeed, assuming hydrostatic equilibrium one can use 
the X-ray image of a cluster to reconstruct from its temperature and luminosity distribution the 
baryonic and dark matter radial profiles, as explained for instance in |52t I53j. The combination of 
all these techniques should allow a robust reconstruction of the dark matter and baryon components 
of the clusters. 

These standard techniques assume that there is no extra contribution to the cluster mass. How- 
ever, the crucial property of c s = models of quintessence is that a third kind of mass, associated to 
quintessence, enters in the game. How are the different tecniques sensitive to these new mass? Lens- 
ing and galaxy velocity measurements will still trace the total mass, now including, after dark matter 
and baryons, also quintessence. On the other hand, optical light associated to galaxies should be 
rather independent of the accreted quintessence. X-ray experiments should indicate a higher gravi- 
tational mass - both from studying the hydrostatic equilibrium and from temperature measurements 
which are sensitive to the gravitational potential well - but the same baryon component. This effect 
is similar to a reduction of the f gas parameter, the ratio between baryonic and dark matter mass. 
Similar considerations apply also to SZ measurements. 

A useful lever arm for distinguishing the accretion of dark energy from the other uncertainties 
in the description of a cluster is the strong red-shift dependence. Quintessence should be relevant 
only at very low redshift. It would be useful to find out the best combination of observables which 
is able to constrain the presence of extra mass in the cluster. Besides a smoking gun of clustering 
quintessence, it would be a useful consistency check between the various measurements. 

8 Conclusions 

Using the spherical collapse model, we have studied how the clustering of quintessence with negligible 
speed of sound can affect the prediction for the mass function of dark halos. As quintessence does 
not develop pressure gradients, it follows geodesic motion remaining comoving with the dark matter. 
In contrast to the case where quintessence remains smooth, spherical collapsing regions behave as 
exact closed FLRW solutions with the quintessence contributing to the overdensity. To study the 
spherical collapse we found it useful to use Fermi coordinates where the effect of the expansion can 
be treated as a perturbation around Minkowski spacetime locally around a spatial point. In contrast 
to comoving coordinates, a unique coordinate system can be employed to describe both the interior 
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and exterior of the collapsing region. 

Quintessence with zero speed of sound modifies dark matter clustering with respect to the smooth 
quintessence case through the linear growth function and the linear threshold for collapse. Besides 
these conventional effects there is a more important and qualitatively new phenomenon: quintessence 
mass adds to the one of dark matter, contributing to the halo mass by a fraction of order (1 + 
w)0,Q/Cl m . This effect is quite difficult to model and in this paper we have adopted a simplified 
treatment which gives an accurate estimate of the high mass tail of the distribution where the effect 
is more relevant. 

As dark energy plays an active role in the formation of structures, the distinction between what 
we call dark matter and dark energy becomes fuzzy. The distinction between the two components 
can be probed by the different redshift dependence. It is quite remarkable that our knowledge of 
structure formation still allows a completely inhomogeneous dark energy component on short scales. 

In this paper we did not attempt to study experimental constraints and forecasts. However, we 
have seen that for values of w which are not too close to the cosmological constant one, let's say 
|1 + w\ > 0.1, the predictions for the high mass tail of the mass function in the c s = and the 
c s = 1 cases are quite distinctive. The effect of clustering quintessence gives order one changes in 
the expected number of clusters. Whether the effect will be measurable or not will depend on the 
possibility of breaking the degeneracy with cosmological parameters, most notably as and f2 m , and 
on the good recostruction of the limiting mass of the survey. Similar concerns have been addressed 
in the context of the effects of primordial non-Gaussianity on the mass function in |54j . 

Our work can be continued in various directions. A better theoretical treatment of quintessence 
accretion would strengthen our predictions for the mass function, although at a certain point only 
a (challenging) numerical simulation can make a fully reliable prediction. On the data side it 
would be interesting to understand what is the minimum value of |1 + w\ that allows a distinction 
between c s = and c s = 1 using the forthcoming cluster mass function measurements. One should 
also explore whether other probes, lensing for example, are better suited to study the non-linear 
clustering of quintessence. It would also be interesting to explore the effect of quintessence on the 
halo dynamics, but this probably requires a way to deal with the formation of caustics. In a more 
model-independent way, one could try to use data on clusters to constrain the presence of extra 
gravitational mass besides the dark matter mass. We leave all this for future work. 
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